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Abstract. - Polarons, introduced by Davydov to explain energy transport in a-helices, correspond 
to electrons localised on a few lattice sites because of their interaction with phonons. While the 
static polaron field configurations have been extensively studied, their displacement is more difficult 
to explain. In this paper we show that, when the next to nearest neighbour interactions are included, 
for physical values of the parameters, polarons can spontaneously move, at T = 0, on bent chains 
that exhibit a positive gradient in their curvature. At room temperature polarons perform a random 
walk but a curvature gradient can induce a non-zero average speed similar to the one observed at 
zero temperature. We also show that at zero temperature a polaron bounces on sharply kinked 
junctions. We interpret these results in light of the energy transport by transmembrane proteins. 



Introduction. — Proteins, essential components of all biological cells are central to 
their proper functioning. As "form determines function" , a study of protein structure and 
dynamics is of utmost importance in elucidating their role in cellular behaviour. One of 
the key problems in biology is to understand energy transport from one part of the cell to 
another and to study the role of protein conformations and conformational transitions in 
this process. 

The mechanism of charge and energy transport in proteins and other bio-macromolecules 
at the atomic scale was proposed Davydov and co-workers [T| . In this approach the transport 
properties are considered in terms of the emergence of a 'polaron' whose properties and 
dynamics are used to describe the resultant transport. The polaron describes a localised 
excitation which carries energy corresponding to some vibrational modes of a group of 
molecules and a distortion of the chain containing these molecules. The system acts as 
a particle and its dynamical properties can be studied in terms of solutions of particular 
differential equations that describe its behaviour. 

The Davydov theory hinges on the assumption that an extra electron or energy quanta 
released in the hydrolysis of ATP (adenosine triphosphate) can be stored by the protein 
molecule in its vibrational mode. The non-linear coupling between the vibrational mode 
and other excitations on the chain leads to the formation of a soliton (polaron in the case of 
the interaction between the phonon modes and the electron in a polarizablc medium). The 
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soliton (polaron in this case) can then propagate along the polypeptide backbone leading 
to energy /charge transport. The soliton mediated transport mechanism has been applied 
to helical proteins [2], The theory of non- linear energy transport in bio- molecules is re- 
viewed in 3,4 . Most of the studies are based on simple models of one-dimensional reduc- 
tions of the three-dimensional structure of the a- helix proteins i.e., only a single strand of 
hydrogen-bonded peptide units is considered. Moreover, the original set of discrete lattice 
equations is often treated in a continuum approximation. Recent studies based on two and 
three-dimensional models describing the solitonic and/or polaronic transport of energy are 
described in detail in [SHE]. 

Studies of polaron transport on a helical polymers have recently been reported by Pfen- 
ning [9] . The Henning model is an extension of semiclassical treatments of polaron transport 
within the Holstein model which incorporates three strands that comprise the a helix. Nu- 
merical schemes and variational approaches have been utilized in the treatments of the 
Holstein model with a hard non-linear on-site potential term in one dimension to obtain 
polaronic ground states, their phase diagram in terms of the non-linearity parameter as 
well as the one characterising the electron-phonon coupling [TUHT2"] . It has also been shown 
that the presence of the non-linear on-site potential polaron formation is possible in higher 
dimensions, a fact borne out by the Henning model. Normal mode analysis of this model 
revealed the existence of a low-frequency pinning mode, and hence enabled the construction 
of one and two-dimensional moving polaron solutions by a Floquet analysis |10[lll|. The 
effects of the lattice non-linearity were reported to have led to a dramatic reduction of the 
polaron effective mass. 

However, in most of these studies the polaron is "kicked" from its rest state via a per- 
turbation |13| . In this paper we look at the spontaneous polaron transport via charge- 
conformational coupling. In particular, we study polaron transport on a flexible chain with 
an imposed initial bend at both zero (T = 0) and non-zero (T ^ 0) temperatures. We now 
summarise our main results. 

At T = a polaron can undergo spontaneous motion via the coupling of charge with the 
conformational degrees of freedom. Thus an imposed bend on the chain causes the polaron 
to accelerate. However, we have found that when such an accelerating polaron encounters 
a kink, a slope discontinuity along the chain backbone, it gets reflected instead of contin- 
uing along its original direction of motion guided by inertia. The reason for this due to 
the long-range nature of the exchange coupling term. At finite temperature T ^ thermal 
fluctuations wash away directed transport and the polaron undergoes large amplitude fluc- 
tuations about its mean position. However, since the chain ends act as reflecting walls, a 
polaron formed closer to one end of a straight chain would often get reflected, resulting in 
an overall small but non- negligible drift . 

The paper is organised as follows: In the next section we present the Hamiltonian describ- 
ing polaron transport on a flexible chain and change variables to dimensionlcss quantities. 
The dimensionlcss parameter values corresponding to physically relevant systems e.g. pro- 
teins are computed in Sec. 3. We describe the results of our numerical studies in Sec. 4 and 
we set our work in perspective in the final section. 

Polaron Hamiltonian on a Flexible Chain. — We consider a modified version of 
the polaron model proposed by Mingaleev et al. |14j . Our model involves a semi-classical 
treatment of the interaction between a phonon field R n and an electron field </>„ on a flexible 
linear chain whose nodes are labelled by an index n. 

The Hamiltonian of the model is given by 

n = E 

n 

where M is the mass of the chain node, W is the linear excitation transfer energy and A 
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the non-linear self-trapping interaction. The excitation transfer coefficients J n _ r 
form: 

Jn,m = J(\fn ~ ? m \) = (e a - 1) e^""^/ , 



are of the 



(2) 



where a -1 sets the relative length scale over which the interaction decreases, in units of 
a, and where a is the rest distance between two adjacent sites. The J n _ m describes the 
long range interaction between the electron field at different lattice sites n and m; its value 
decreases exponentially with the distance between them. 

Note that the normalisation of the electron field is preserved in our model, i.e. 
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The phonon potential U n consists of three terms: 
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The first term in Eq.Q models the elastic energy describing the stretching of the adjacent 
nodes of the chain, where a is the equilibrium separation between them. The second term 
describes the bending energy of the chain akin to semi-flexible polymers. Here Lp n is the 
rest angle between adjacent lattice links and 9 max is the largest angle allowed. In the 
Mingaleev model [14] ip n = and the equilibrium configuration of the chain is a straight 
one. Non-zero values of tp n lead to a bent chain. The last term, (proportional to 5), models 
hard-core repulsion between the atoms of the chain. However, we note that this term does 
not contribute significantly in determining the equilibrium chain conformation. 

In this paper symbols denoted by an overhead carat sign e.g. M, a etc. correspond to 
physical variables carrying units and dimensions while those without it correspond to non- 
dimensional variables and parameters. We rescale time r by defining a timescale tq = A^s 
and rescale distances by a length scale a. The non-linear coupling parameter g appears as 
a dimensionless ratio of two energy scales. 
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In terms of these variables the Hamiltonian takes the form 
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Writing r n — (xi jTl , X2. n , X3,m ) we can derive the equation of motion for x^ n from the 
Hamiltonian Eq.([B]). Thermal fluctuations are incorporated by adding a delta correlated 
white noise F(t) which satisfies 

< F(0)F(s) >= 2Tk B TS(s). (9) 

The resulting Langevin equations describing the coupled dynamics of the chain and of the 
polaron is given by 

m fe rn<k *' 

i ~tt L _ + ^ Jnm< ^ m + f l^nl 2 ^ = 0, (10) 

Note that we also have 

W 2 

k B T = k B T— = k B TWg (11) 

and, as the equation for Xi is expressed in units of W 2 /(Aa), we have T = Ta 2 /H. 

Physical Parameter Values. — In order to study the feasibility of the mechanism 
of the energy transport via polarons for biologically relevant systems we have used the 
parameter values corresponding to those of a-helices. 

We note that for Amid-I vibrations in a-helices [TS] we have: W ~ 2 x 1Q~ 22 J w l.2meV, 
a = 19.5N/m and M = 2 x \Q- 25 kg. a sa 0.45nm. e = 0.02eK A = 8^ = 4.74 x 1CT 22 J = 
O.OOSeV^. Moreover, k can be evaluated from the persistence length of a-helices A « 65nm 

k = \k B T/a « 6 x 10 -19 J. (12) 

Though we do not have experimental values of <5, it is clear that it must be larger than a. 
For the friction coefficient we assume that T m 6ir fj,a, where /i = 0.001 Pas is the water 
viscosity. Incidentally, for the cytoplasm \i is up to 4 times larger than this value. 
At "room" temperature T = 300-fT, the non-dimensional parameter values are 

.9 = 2.5 o = 51351. fc = 7776 T = 8143 (13) 

M = 2.79 x 10 5 k B T = 54 v = 0.01 r = 10" 12 s. 

Following Mingalccv et al. [14] we choose a = 2. 

Next we have performed the simulation of the time evolution described by the equations 
Eq. (|10[) - to explore the coupled dynamics of polaron transport and the conformational 
transitions of the chain associated with such motion. For the values of g and k mentioned 
above, we have found that the physical a-helices are too rigid for the polaron to bend the 
chain and move along it as suggested for DNA in [T3] (see Fig. 2 in [T3] where our g 
corresponds to N). 

As a matter of fact, physical polarons always have a relatively small energy, too small 
to bend a polymer spontaneously, even for DNA. Nevertheless, the Mingaleev et al. model 
can be used to study the properties of polarons on chains that arc, like most proteins, 
naturally bent. One can expect the electron-phonon interaction to lead to the spontaneous 
displacement of polarons on bending gradients for the following reason: the main effect of 
the electron-phonon interaction is to favour configurations where the distance between nodes 
is reduced (i.e. the interatomic separation is smaller than the equilibrium lattice spacing 
in the vicinity of the polaron). As a kink in the linear chain reduces the distance between 
nodes that are not adjacent to each other, one expects that the polaron would be attracted 
by the kink. Thus chain conformations where the angle between consecutive tangent vectors 
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increases monotonically serve as an attractive potential as seen by the polaron. It is thus 
expected that the polaron will accelerate towards a region of high curvature. 

In the next section we present the results of investigations of whether these expectations 
are correct for the physically relevant proteins (i.e for chains with physically relevant values 
of their parameters). 

Polaron on a Bent Chain. — We consider a chain of N = 60 points with a bent 
mid-section i.e. the region between n — 25 and n = 45. Such a configuration is achieved by 
setting, initially, <p n as follows: 

n < 25 tp n = 

25 < n < 45 ip n = (n - 25)dtp (14) 
n > 45 ip n = 0, 

where dip is the incremental increase in the angle between adjacent links. 

In order to numerically simulate the dynamics of the polaron on such a chain, we initially 
generated and saved the polaron on an undeformed (i.e. straight) lattice (ip n = for all 
n). Such a polaron was obtained by a relaxation method: a friction term was added to the 
phonon field and the eigenvalue problem for the electron field was solved simultaneously by 
relaxation. Starting with an electron spread over a few lattice sites and located on a node 
which we have chosen to be n = 30, the system was relaxed until both the electron and the 
phonon fields reached a stationary configuration. This produced a polaron configuration 
centred at n = 30. The lattice was then deformed by changing the values of p n to those 
of (|14p and a second relaxation was performed to let the phonon field reach its equilibrium 
configuration. The electron field of the relaxed polaron was then restored, the time set to 
t = 0, and the equations for the phonon and electron fields integrated numerically. 

The zero temperature (T = 0) dynamics was obtained by numerically integrating the 
electron and phonon fields in Eq. (|10l) employing a fourth order Rungc-Kutta scheme, noting 
the time it took for the electron to start moving. 

Fig. [^summarises our results for T — 0. Panel (a) shows typical trajectories of a polaron 
for different values of a, the parameter that models the range of the long-ranged interaction 
Jmn (Eq.©) Our simulations show that the main effect of including the next to the nearest 
neighbour interaction is to make the polaron "to get attracted" to the bend. This is caused 
by the fact that the lattice points near this bend are nearer to each other and this lowers the 
energy of the polaron. Hence, on a bent chain, with a slowly increasing bending angle, the 
polaron is attracted increasingly to the bend as it moves along the chain. This corresponds 
to the spontaneous displacement of the polaron. It is curious to note that upon reaching 
the edge of the bent chain beyond which the chain is straight the polaron is reflected. This 
can be explained by the fact that the energy of the polaron is the lowest where the bending 
angle is large and the highest where the chain is straight. The transition from a bent to a 
straight configuration thus corresponds to a potential wall on which the polaron bounces. 
One notices also that as the parameter a is increased the exchange energy between non 
adjacent nodes decreases more rapidly with the distance separating them. This results in a 
smaller acceleration of the polaron as seen in Fig. 1. The speed of the polaron was evaluated 
as V = [n(to + At) — n(to)]/At where to is the time at which the polaron started to move 
and At = 50. When the polaron travelled a distance larger than 20, we took V = 20/ At 
where At was the time taken to cover this distance. 

Fig. [Ub) shows the chain configuration with the electron density superimposed on it. 
Each filled circle (black) corresponds to a lattice point, while the electron density, denoted 
in grey /yellow, is overlayed on them. The intensity of the lighter circles provides a measure 
of the electron density at a particular site, i.e. a lighter colour corresponds to a higher 
electron density. The maximum in the electron density corresponds to the position of the 
polaron used to generate Fig. [2(a). 
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Fig. 1: Polaron displacement on a bent chain with dip — 0.01 at T = for a TV = 60 chain with 
the initial polaron position at xo — 30. The time step of integration is chosen to be dt = 50. Panel 
(a) shows the trajectory of the polaron on the bent chain for different values of a (a = 1 (red solid 
line), q = 1.5 (blue dashed line), a = 2 (green dash-dotted line), and a — 2.2 (magenta dotted 
line). Panel (b) shows the chain configuration and the electron density. Each circle corresponds to 
a node and a lighter colour corresponds to a higher electron density. 



Next we have investigated the dependence of the polaron dynamics on dip - the gradient 
of the bending angle of the chain. Fig. [2] shows the variation of the average velocity (V) 
as a function of dip. For our choice of parameter values the polaron does not move until a 
critical value dip m 0.0056 is reached. The dependence of the velocity on dip in the critical 
region follows a power law (V) ~ (dip — dip c )" akin to elastic depinning of interfaces. In 
order to express the average speed in physical units, we have to multiply the dimensionless 
speed by a/t « 328.5^. Thus the average speed of the polaron (V) « IQnm/ns. 

In order to study the effect of thermal fluctuations on the polaron transport we have 
performed finite temperature T / simulations. For such simulations we generated a 
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Fig. 2: Figure showing the variation of the velocity (V) as a function of the bend angle dip for a 
bent chain of TV = 60 nodes. The initial position of the polaron is at xo = 30, the bent region 
extends from n = 25 to n = 45 and T = 0. The speed was computed as V — d/At where At = 50 
and d is the distance travelled by the polaron during this time interval. 
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Fig. 3: Figure shows the variation of the average velocity (V) as a function of the bend angle dtp 
for a N — 60 chain. The chain is bent between n — [25, 45] nodes, and the initial position of the 
polaron is at xq — 30. The temperature T = 300K, and the data is averaged over N run = 1000 
simulations. The speed was evaluated exactly like in figure 2. Inset shows the variance of the 
velocity vs. dip data. 

polaron on a straight lattice at T = as above and saved it. Next, we deformed the lattice 
into a bent configuration as above and integrated our equations for 300 units of time to ensure 
thermal equilibrium. We then restored the electron field, without altering the phonon field, 
and numerically integrated Eq (|10j) with the noise term to investigate the motion of the 
polaron along the chain. This simulation strategy, faithfully mimics the sudden excitation 
of a polaron and its time evolution. 

Our finite temperature studies have focused on the polaron dynamics at "room tem- 
perature", i.e. T = 300-ftT. Clearly, the effects of taking non-zero temperature are very 
significant as seen in Fig. [3] It is clear from this figure that the polaron undergoes large 
amplitude fluctuations about its initial position and its velocity is significantly altered by 
the thermal effects. The thermally averaged polaron velocity shown in Fig. [3] is averaged 
over N run = 1000 simulations, each time having measured the average displacement over 
t = 50 units of time. In all cases the life time of the polaron was of the order t = 200, in 
our units, which corresponds to about 300ps. Afterwards the electron was dclocalised on 
the lattice. In Fig. |4] we have plotted the variation of the thermally averaged velocity as a 
function of dtp; this is to be compared with the curve obtained for T = in Fig. [5J 

Our simulations clearly show that the thermal effects lead to a wider displacement of 
the polaron wich can involve motion in both directions. The polaron is thus subjected to a 
random walk motion in a system with an attractive force provided by the bend of the chain. 
So, while the displacement of a polaron in a protein at room temperature is random, its 
position is still biased by the bending gradient; a thermalised polaron thus follows a random 
walk, but its average displacement is similar to that of a polaron at T = 0. 

In order to quantify the dependence of the bend on the spontaneous transport of polarons 
we have calculated the dependence of the critical incremental angle between segments dip c 
on its initial distance from the bend at T = 0. This is shown in Fig. |4] It is clear from the 
plot that dip c oc exp[— where £ is a constant and L is the distance between the polaron 
and the edge of the region with the bend. 

Finally we have also explored a possible mechanism of the energy /polaron transport for 
biomolecules i.e. a helices in our case. Thus we have considered the transport of polarons 
on straight chains at finite temperatures. Since the edges of the chains act as reflecting walls 
a polaron that is initially formed close to one edge of the chain will get reflected from it 
during the course of its random walk motion with an amplitude set by the thermal energy 
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Fig. 4: Figure shows the variation of the critical bend angle dtp c as a function of the initial relative 
polaron position L with respect to the bent region, no lPO ; = 25 + L. N = 60 and T = 0. The speed 
was evaluated exactly like in figure 2. 



scale. Fig. [5] shows the dependence of the mean position of such a polaron that has initially 
been created at xo = 5, 10, . . .25. While the steady state polaron position is the same for 
all values of xq, their early time behaviour is significantly different. Thus a polaron formed 
at xo = 5 has a high average velocity due to its reflections from the n = edge and is, 
on average, transported further at a given time (say t = 140 ) than a polaron formed at 
x a = 25. This compels us to conjecture the following model of polaronic transport for 
biophysical systems: if a polaron is generated by electron-phonon interactions near an edge 
of a short straight segment of a protein it will undergo reflections from the nearby edge, 
leading to the directed motion towards the other edge. When the system contains receptor 
molecules that can absorb this non-linear excitation and transfer it to other molecules or 
polymers, the polaron displacement can be used to transport the energy released by ATP 
hydrolysis to a nearby regions of the cell where it is needed. 
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Fig. 5: Figure shows the average position of a polaron as a function of time for a straight chain 
of N = 60 segments at T = 300A' for different initial positions (xi = 5 (red diamonds), x% = 10 
(black circles), x% = 15 (blue squares), x% = 20 (magenta upper triangles) and xt = 25 (green stars) 
) averaged over N run = 1000. Inset shows the average position and the fluctuation about the mean 
position of the polaron for two initial positions Xi — 5 and Xi = 25. 
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Conclusions. — In this paper we have studied the transport of energy along bio- 
polymers via polaronic mode. Starting from a model that describes coupled electron phonon 
dynamics, we chose physically relevant parameter values and studied polaron dynamics on 
bent configurations of bio-polymers. We have found that the bent chain induces spontaneous 
polaron transport due to some of its nodes being closer together and so generating an 
attractive potential for the polaron. When we looked at this problem for proteins at zero 
temperature the movement of polarons was observed when the gradient of the bend was quite 
small (dip < 0.0055 in our units). The polaron moved at a speed V « Vonmjns and lived for 
about t = 300ps. The thermalisation of the polaron, at room temperature, modified some 
aspects of the polaron transport: the polaron now exhibited a random walk motion biased 
in the direction of the bending gradient. The bending gradient of the chain thus induced an 
average polaron displacement. Moreover, on a straight chain, we observed that a polaron 
near the edge of the chain exhibits spontaneous displacement in the direction away from 
the edge. We have thus shown that polaron displacement can be induced spontaneously 
without the need to kick the polaron in one direction and that the direction of transport 
can be determined by the configuration of the polymer. 
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